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ABSTRACT 

We present a newly observed relation between galaxy mass and radial metallicity gradients of early- 
type galaxies. Our sample of 51 early- type galaxies encompasses a comprehensive mass range from 
dwarf to brightest cluster galaxies. The metallicity gradients are measured out to one effective radius 
by comparing nearly all of the Lick absorption- line indices to recent models of single stellar populations. 
The relation shows very different behavior at low and high masses, with a sharp transition being seen 
at a mass of ~ 3.5 x 10 10 M Q (velocity dispersion of ~ 140 km s _1 , Mb ~ —19). Low-mass galaxies 
form a tight relation with mass, such that metallicity gradients become shallower with decreasing 
mass and positive at the very low-mass end. Above the mass transition point several massive galaxies 
have steeper gradients, but a clear downturn is visible marked by a broad scatter. The results are 
interpreted in comparison with competing model predictions. We find that an early star-forming 
collapse could have acted as the main mechanism for the formation of low-mass galaxies, with star 
formation efficiency increasing with galactic mass. The high-mass downturn could be a consequence of 
merging and the observed larger scatter a natural result of different merger properties. These results 
suggest that galaxies above the mass threshold of ~ 3.5 x 10 10 M Q might have formed initially by 
mergers of gas-rich disk galaxies and then subsequently evolved via dry merger events. The varying 
efficiency of the dissipative merger-induced starburst and feedback processes have shaped the radial 
metallicity gradients in these high-mass systems. 

Subject headings: galaxies: dwarf — galaxies: elliptical and lenticular, cD — galaxies: formation — 
galaxies: evolution — galaxies: stellar content 



1. INTRODUCTION 

Radial metallicity gradients have been found to charac- 
terize the stellar population of early-type galaxies. They 
are the chemodynamical fossil imprints of galaxy forma- 
tion and evolution mechanisms. The presence of such 
population gradients and their relationship to galaxy 
structural parameters provides a strong constraint on 
galaxy formation. Competing scenarios that describe the 
formation of early-type galaxies predict contrasting be- 
haviors for the mass-metallicity gradient relation. 

In the classical models of monolithic collapse, stars 
form during a rapid dissipative collapse while the gas 
sinks to the c enter of the forming galaxy (lLarsoiil 1 1974t 
iLarsonl Il97l : iCarlberd fl98i lArimoto fc Yoshiil 11987ft . 
The inflowing gas is chemically enriched by evolving 
stars, and contributes metal-rich fuel for star formation. 
Thus a negative radial gradient is established, in which 
central stars are more metal-rich than those born in outer 
galaxy regions. Metallicity-dependent gas cooling and a 
time delay in the occurrence of galactic winds (i.e. on- 
sets that vary with the local escape vel ocity) cooperate 
in steepening any me t allicity gradients ( Matteuccllll994 
iMartinelli et al.lll998t[Pipino et al.ll2008l) . TJie^redicted 
metallicity gradients are as steep as -0.35 ()Larsonlll974l ). 

1 Based on observations obtained at the Gemini observatory, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under a cooperative agreement with the NSF 
in behalf of the Gemini partnership: the National Science Foun- 
dation (United States), the Particle Physics and Astronomy Re- 
search Council (United Kingdom), the National Research Coun- 
cil (Canada), CONIC YT (Chile), the Australian Research Council 
(Australia), CNPq (Brazil), and CONICET (Argentina). 

2 Corresponding Author: mspolaor@astro.swin.edu.au 



-1.0 (lLarsonlll975h . and -0.5 dex per radius dex ([Carlberd 
Il984h . The depth of the galaxy potential well regulates 
the effectiveness of these processes, so that a strong cor- 
relation between metallicity gradient and galaxy mass 
is expected. High-mass galaxies retain more metals al- 
lowing steep negative metallicity gradients of -0.5 dex 
per radius dex, co mpared to low-mass systems with al- 
most zero gradient (lGibsoni ri997: Chi osi fc Carrarol2002l : 
IKawata fc Gibson! I2003D 

In the hierarchical clustering scenario of galaxy for- 
mation, early-type galaxi e s are t he product of merging 
event s jK auffman n et all 119931: iKa uffmann & C harlotl 
[19981: IDe Lucia et all 12001 ; iDe Lucia fc Blaizotl [2001 . 
The turbulent mixing caused by merging events is ex- 
pected to "wash out" any possibl e ambient g radient 
presents in the progenitor galaxies (IWhitd[l980l ). Sim- 
ulations suggest that stellar population and kinematic 
properties of a merger remnant are determined by the 
progenitors' mass ratio an d by the amount o f gas involved 
in the merger event (e.g. JNaab et al.lf2006l ). Metallicity 
gradients are produced by a rapid central merger-induced 
starburst and their strength is p roportional to the effi - 
ciency of the dissipative process (jHopkins et al.ll2008af ). 
The gradients are predi cted to only weakly depend on the 
remnant galactic ma ss ( Bck ki fc Shio^[l999T iKobavashil 
[200llHoT)kins et al.ll2008bl ). 

Metallicity gradients in early-type galaxies have been 
studied by means of increasing strength of metal- 
sensitive lines and the reddening of colors. Unfortu- 
nately, most estimates are made under the crucial a priori 
assumption that this increase is entirely caused by a ra- 
dial variation of stellar metallicity and therefore neglect- 
ing age effects. Typically, past works have derived total 
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metallicities using just a few metal-lines (i.e. Hj3, Mg2, 
^5270) and have sampled the relatively central galactic 
region inside 0.5 effective radii. To date, a clear and con- 
sistent picture of the behavior of the mass-metallicity 
gradient relation has not emerged from observational 
studies. 

In a sample of 42 early- type galaxies, ICarollo et al.l 
(1993) observed that their Mg% gradients became shal- 
lower for galaxies with masses less than about 10 11 M Q , 
while more massive syste ms show no correlation . From 
a literature compilation, iKobavashi fc Arimotol (|1999f ) 
studied metal-line index gradients of 80 elliptical galax- 
ies without findi ng any statistic a lly sig nifica nt trend with 
mass. Recently, lOgando et al.l (|2005l ) and iForbes et al.1 
(2005) claimed opposite trends fo r a sample of 3 5 and 
29 massive early-type galaxies. lOgando et al.l (|2005l ) 
showed that both flattening and steepening of metallicity 
gradie nts are possible in m assive galaxies. On the other 
hand, IForbes et al.1 (|2005f ) found that gradients become 
steeper in more massive systems. 

Here we present new results from high signal-to-noise 
Gemini GMOS observations of 14 low-mass early-type 
galaxies in the Virgo and Fornax clusters (M. Spolaor 
et al., in preparation). By u sing all of the Lick indice s 
and adopting the technique of lProctor fc Sansoml ([2002) , 
we have been able to break the age-metallicity degen- 
eracy and obtain reliable metallicity gradient measure- 
ments out to more than two effective radii. Combin- 
ing our measurements with those from a novel literature 
compilation, we extend our study to a broad mass range 
(i.e., from dwarfs to brightest cluster galaxies). The final 
sample of metallicity gradients in 51 early- type galaxies 
is a uniquely uniform sample due to the use of nearly all 
of the Lick indices and the steps taken to separate the 
age-metallicity degeneracy. 

In this Letter, we are able to probe for the first time 
the low-mass end of the mass-metallicity gradient rela- 
tion. A clear picture of the relation along the full mass 
range is finally possible. Its interpretation in the con- 
text of competing model predictions allows us to place 
constraints on galaxy formation mechanisms. 

2. A CONSISTENT DATA SAMPLE 

Long-slit spectroscopic observations of 14 dwarf Virgo 
and Fornax cluster early-type galaxies were performed 
with GMOS at the Gemini South telescope during 
the semesters 2006B (Program ID GS-2006B-Q-74) and 
2008A (Program ID GS-2008A-Q-3). The galaxies were 
selected to be low-mass using the central stellar velocity 
dispersion (a) and the total B-band absolute magnitude 
(Mb) as independent proxies of galaxy mass. The sam- 
ple uniformly covers the range 1.6 < loger < 2.15 and 
— 19.5 < Mb < —16. The spectrograph slit was set to 
0.5 arcsec and covered 5.5 arcmin in length. We used 
the 600 line mm -1 grism blazed at 5250 A providing 
a dispersion of 0.458 A per pixel and a spatial resolu- 
tion of 0.1454 arcsec per pixel. The spectroscopic data 
were reduced using the Gemini data processing software. 
We spatially resolved several apertures out to one effec- 
tive radius along the major axis of each galaxy. The 
spatial width (i.e. the number of CCD rows binned) for 
each extracted aperture increases with radius to achieve 
a signal-to-noise ratio of ~ 30 at the central wavelength 
of 5200 A. The spectra cover the full range of the 21 line- 



strength indices of th e original Lick/ IDS (Image Dissec- 
tor Scanner) system (jWorthevlll994l ) plus the additional 
age-sensitive Balmer indices HSa, H^a, H8f and H^p 
(|Worthev fc OttavTa ni 1997). The measured indices were 
calib rated to the L i ck syst em as described in Appendix 
A of ISpolaor et al.1 (|2008b[ ). We note that a comprehen- 
sive analysis of line-strength indices of these new data 
are the subject of a future work by M. Spolaor et al. (in 
preparation) . 

Independent radial profiles of luminosity-weighted age, 
total metallicity [Z/H] and a-abundance ratio [a/Fe] 
were derived fro m the fully calibrat e d indi ces. Briefly, 
the technique of iProctor fc Sansoml (|2002f ) is a statis- 
tical x 2 minimization of the deviation between the ob- 
served and modeled index values as a fraction of the in- 
dex errors. The strength of the method is that it works 
with as many indices as possible in order to break the 
age-metallicity degeneracy that affects each index differ- 
ently. Consistent with the other studies, we u s ed th e 
single stellar population models of iThomas et al.l (2004), 
which include the variation of Lick indices to abundance 
ratios. 

Gradients were obtained by a linear least-squares fit to 
the radial metallicity profile weighted by the errors asso- 
ciated with each metallicity value. Therefore, metallic- 
ity gradients express the derivative variation of the total 
metallicity along the galactocentric radius (scaled by the 
effective radius). We represent the values in logarithmic 
space as A[Z/H]/A\og(R/R e ) (hereafter A[Z/H}). To 
be able to compare gradients from different data samples 
we uniformly fit the data points beyond the seeing limit 
(i.e., > 1 arcsec) out to one effective radius. We note 
that the inclusion of these very inner regions introduces 
only a statistically insignificant change in the gradients. 

Our full sample is a compilation of 51 metallicity gra- 
dients in early-type galaxies from six different studies 
that adopted the same criteria and analysis technique. 
Spolaor et al. (2008a, 2008b) studied the star formation 
and chemical evolutionary hist ory of the two brightest 
galaxi es of the NGC 1407 group. ISanchez-Blazquez et al.l 
(|2007| ) used Keck LRIS long-slit spectra to study a 
sample of 11 early- type galaxies covering a wide range 
in mass, situ ated in the field, p oor groups and the 
Virgo cluster. iBrough et al.1 (|2007f ) focused on the stel- 
lar population and kinematic properties of three nearby 
brightest group galaxies a nd 3 brighte s t clus ter galax- 
ies located at z ~ 0.055. iReda et ah! (2007) analyzed 
the metallicity gradients of 12 isolated luminous early- 
type galaxies. The sample of Proctor (2003) consists 
of 11 early- type galaxies located in the Leo cloud and 
the Virgo cluster. We note that five of the 11 galax- 
ies of the ISanchez-Blazquez et al.l (|2007t ) sample over- 
lap with the Proctor (2003) sample; their gradients are 
similar, however we adopted the Proctor values due to 
his use of all 25 Lick indice s versus the 20 used by 
ISanchez-Blazquez et all (|2007f ). 

3. RESULTS AND DISCUSSION 

We measured the central stellar velocity dispersion for 
our galaxy sample by averaging the dispersion values 
within a radius of R e /8. Total B-band absolute mag- 
nitudes were estimated from total apparent magnitudes 
given in the NED database and corrected for Galactic ex- 
tinction. For each galaxy we assumed the distance mod- 
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Fig. 1. — Metallicity gradients A[Z/H] as a function of galaxy central velocity dispersion logo - (left panel) and total B-band absolute 
magnitude Ms (right panel), both independent proxies of galaxy mass. In the left panel, the 51 data points are coded by original data 
sample. In the right panel, the points are coded by environment: isolated/field, grou p, cluster and brightest group/cluster galaxies (BGG, 
BCG). Also shown are the m ass dependent prediction s from the merging models of Bekki & Shioya (1999) as dot-dashed lines, and the 
dissipative collapse models of Kawata & Gibson (2003) as dashed lines. The region occupied by the remnants of major mergers between 
gas-rich disk galaxies, as simulated by Hopkins et al. (2008a), is shown by the yellow shading. 



ulus f rom s urface brightness fluctuations of iTonrv et alJ 
(120011) andlMei et al.H2005h . We applied the correction 
ofO ensen et al.l (|200*3T ) to the values of Tonry. 

In Figure 1 we present the relation between metallicity 
gradients and two independent proxies of galaxy mass. 
Gradient values of the 51 galaxies are shown as a function 
of (j (left panel), coded by the data source, and Mb (right 
panel), coded by environment. We defined four environ- 
mental classes: isolated/field, group, cluster and bright- 
est group/cluster galaxies (BGGs, BCGs). The results 
are compared to theoretical predictions of the com peting 
galaxy formation mod els oflBekki fe Shioval (11999ft (dot- 
dashed black line) and lKawata fc Gibson! (|2003( ) (dashed 
black li ne). For compl e teness, we also show the simula- 
tions of I Hopkins et alJ (|2008af ) of gas-rich major merger 
remnants which have "excess" central light due to a 
merger- induced central star burst. The region they oc- 
cupy is shown by the yellow shadi ng. 

The models of iBekki fc Shioval (|1999fl simulated the 
formation of early-type galaxies by merging gas-rich disk 
galaxies of varying mass ratio at redshift z > 2. The 
chemodynamical evolution of mass-dependent galaxy 
properties was investigated by inclusion of supernova 
feedback and time delay between star formation and 
metal ejection. As can be seen in Figure 1, the metallic- 
ity gradients found in merger remnants are shallow and 
only weakly de pendent on the remnan t galaxy mass. 

Nevertheless, iHopkins et al.l (|2008aD have been able to 
reproduce merger remnants with a variety of metallic- 
ity gradients (0.0 < A[Z/H] < —0.6) in simulations 
of gas-rich disk galaxies major mergers, including feed- 
back processes from g alactic winds and black hole growth 
(Springel et al.l l2005f h Similarly, they did not find a 
strong dependence with mass but they stressed that the 
degree of dissipation in the central merger-induced star- 



burst is the most important factor in determining the 
gradient strength in the remnant. The impact on stellar 
population gradients of possible subsequent dry mergers 
of these systems has the same effect observed for pas- 
sively evolved gas-rich merger remnants, i n which metal- 
licity g radients tend to flatten with time (| Hopkins et al.l 
I2008bh . 

In the ch e mod ynamical simulations of 
iKawata fc Gibsonl (|2003f ). an elliptical galaxy be- 
gins via the dissipative collapse of a slow-rotating, 
overdense, top-hat protogalactic gas cloud, on which 
CDM perturbations are superimposed. They find a 
mass-dependence of both the supernova-driven galactic 
winds and the gas infall rate, underlying their important 
role in creating and shaping metallicity gradients. As a 
result, metallicity gradients are steeper in higher mass 
galaxies than those of lower mass systems. 

In Figure 1 our new results populate the previously un- 
explored low-mass branch of the parameter space, while 
the literature values sample the high-mass region. Galax- 
ies with loger < 2.15 (Mb > —19) exhibit a tight trend 
in which the metallicity gradient becomes less negative 
and hence shallower with decreasing mass. At the very 
low-mass end (loger < 1.8, Mb > —17.5) gradient values 
ch ange sign becomin g posit ive or null. The predictions 
of IKawata fc Gibsonl (|2003l ) appear to describe the ob- 
served trend, suggesting an early star-forming collapse 
as the main mechanism for the formation of low-mass 
galaxies. Specifically, star formation efficiency increases 
with galactic mass. 

A significantly different behavior can be seen in galax- 
ies with log (7 > 2.15 (Mb < —19). A downturn from the 
tight low-mass trend is visible marked by a broad scatter 
such that galaxies with increasing mass are characterized 
by shallower metallicity gradients. However, some of the 
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Fig. 2. — Metallicity gradients A[Z/H] as a function of galaxy dy- 
namical mass Md yn . The points are coded by the central age value 
of the galaxies: young (0 — 3 Gyr), intermediate (3 — 8 Gyr) and old 
age (> 8 Gyr). Also shown are the mass-metal licity gradient pre- 
dictio ns from the dissipative collapse models of Kawata & Gibson! 
ll2"003h as dashed lines . The models of IBekki k, Shioval (1199911 and 
IHopkins et al. ( 2008a) are not available. 

most massive galaxies classified as BGGs and BCGs have 
steeper negative gradients. The downturn and the broad 
scatter, could be a consequence of merging and of the 
amount of gas involved in the events, which characterize 
the formation of more massive galaxies. If a significant 
fraction of gas dissipates in a rapid central starburst, 
then, the metallicity gradient of the remnant will reflect 
the mean difference between that of the old pre-merger 
population and that set by the newly formed young stars. 
The gradients of the pre-merger populations have been 
partially flattened by violent relaxation, which acts most 
strongly at larger galactocentric radii. The metallicity 
gradient in the remnant is expected, initially, to be steep 
and then to weaken sli ghtly over ~ 3 Gyr after the merger 
(|Hopkins et al.ll2008al) . 

In Figure 2 we show metallicity gradient values as a 
function of dynamical mass. We estimated the dynamical 
mass of the galaxies by log(M ( j yn ) = 2 log(cr) + log(i? e ) + 
3.1, where R e is the effective radius expressed in pc. By 
averaging the ages estimated from the models within a 
radius of R e /8 we obtained the galaxy central age. In 
Figure 2 we coded the points by assigning three age bins: 
0-3 Gyr, 3-8 Gyr, > 8 Gyr. 

The tight low-mass end of the relation is popu- 
lated by an almost even number of galaxies with 
young/intermediate (0 — 8 Gyr) and old central ages (> 8 
Gyr). The galaxies with younger ages have positive or 
null gradients. Above the high-mass turnover the number 
of galaxies with old central age becomes prominent. We 
quantified the mass transition point around the dynami- 
cal mass of ~ 3.5 x 1O 1O M0, w hich is similar t o the value 
of lO n M found in the data of lCarollo etaD lfl993h . We 
note that the mass value of ~ 3 x 10 10 M© is recurrent 



in literature as critical transition mass for physical prop- 
erties of early-type galaxies, suc h as feedback processes 
from active galactic nuclei (e.g., iKauffmann et all 120031 : 
ICroton et alj|2006t ICattaneo et al.ll2008f) . 

A broader interpretation is that processes of both for- 
mation scenarios arise in the CDM framework for galax- 
ies more massive than ~ 3.5 x 10 10 M Q (loger > 2.15, 
Mb < —19). A rapid star- forming collapse might have 
acted at primordial times creating low-mass progenitor 
galaxies which have then merged to form today's early- 
type galaxies. Therefore, merger events are the main 
mechanism acting to grow galaxy mass and in shap- 
ing the subsequent evolution of galaxies above this mass 
threshold. This has already been suggested by the theo- 
retical work of iKobavas hi (2004). In her chemodynami- 
cal simulations the metallicity gradients do not correlate 
with the galaxy mass but can be used, in principle, to 
infer the galaxy merging history. 

The significant increase in scatter observed beyond the 
high-mass turnover could be caused by the different de- 
gree of dissipation that characterize merger-induced star- 
bursts and therefore the strength of the remnants' metal- 
licity gradients. A further contribution to the scatter 
might come from a time-evolution of the merger-induced 
metallicity gradients. 

4. CONCLUSIONS 

For the first time we have probed the relation between 
galaxy mass and radial metallicity gradients of early- type 
galaxies. Our results show a sharp transition between the 
low and high-mass regimes of the relation at a galaxy 
mass of - 3.5 x 10 10 M Q (log a ~ 2.15, M B ~ -19). 

Galaxies with masses smaller than this transition mass 
are observed to form a tight relation. Their metallic- 
ity gradients become less negative and hence shallower 
with decreasing mass. The interpretation of the results in 
comparison with theoretical models suggests that these 
low-mass galaxies formed during an early star-forming 
collapse, with star formation efficiency increasing with 
galactic mass. 

Above the high-mass turnover of ~ 3.5 x 10 10 M Q 
several massive galaxies have steeper gradients, but a 
clear downturn from the tight low-mass trend is visible 
marked by a broad scatter such that gradients become 
shallower with increasing mass. The downturn could 
be a consequence of merging, and the observed scatter 
a natural result of merger properties. The results sug- 
gest that galaxies formed by gas-rich mergers and then 
evolved, growing their mass via subsequent dry merger 
events. The merger-induced starburst and feedback pro- 
cesses are responsible for creating the radial metallicity 
gradients which are conserved during dry mergers, but 
tend to flatten with time. 
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